General formalism of Hamiltonians for realizing a prescribed evolution of a qubit 
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We investigate the inverse problem concerning the evolution of a qubit system, specifically we 
consider how one can establish the Hamiltonians that account for the evolution of a qubit along a 
prescribed path in the projected Hilbert space. For a given path, there are infinite Hamiltonians 
which can realize the same evolution. A general form of the Hamiltonians is constructed in which one 
may select the desired one for implementing a prescribed evolution. This scheme can be generalized 
to higher dimensional systems. 

PACS numbers: 03.65.Yz, 03.65.Vf 

INTRODUCTION 

Evolution of a quantum system is completely deter- 
mined by the given dynamics. For a closed system, the 
density matrix p{t) of the system satisfies the evolution 
equation, 

U{t)p{Q)U{t)\ 



Pit) 



(1) 



where U{t) — T ex-p[—i J^ H{t)dt\ is the unitary operator 
corresponding to the Hamiltonian H(t). For an open sys- 
tem, its evolution is generally nonunitary and therefore it 
cannot be described by the above equation. However, an 
open system, denoted as Sa, can always be regarded as a 
part of a larger closed system, denoted as Sat, comprising 
Sa and its environment, denoted as St- The evolution of 
combined system Sat now obeys Eq. (^ whereas the evo- 
lution of the open system Sa is described by the Kraus 
representation 1] 



p(i) - ^ A/^p(0)Af^t^ 



(2) 



where Kraus operators Af^t satisfy J2 -^^u-^^/j ~ ^■ 

The Kraus operators are usually expressed as M^ ~ 
{fi\Uab{t)\0), where Uab{t) = T exp[-i J* Hab{t)dt] is de- 
termined by the Hamiltonian Hab{t) of the combined sys- 
tem and l/z) are the orthogonal bases of the environment. 
Clearly, given the Hamiltonian H{t) for the closed sys- 
tem or the Hamiltonian Habit) for the combined system, 
the evolution of the density matrix p{t) is completely de- 
termined through Eq. JQ or Eq. Q. 

It is interesting to study the inverse problem. Sup- 
posing the density matrix p{t) of a quantum system, be 
it closed or open, is given as a time-dependent function, 
how can one obtain the Hamiltonian? This is an inter- 
esting and nontrivial issue because oftentimes in exper- 
iment, one needs to prepare a quantum system which 



is expected to evolve along a given path on the surface 
or inside the Bloch sphere. Certainly, one can always 
find this Hamiltonian through trials and errors. To this 
end, one may look for some trial Hermitian operators and 
require them to satisfy the evolution equation with the 
given p(t), however, this is not easy. Indeed, it may hap- 
pen that even if the form of the Hamiltonian is eventually 
found, it may still be too difficult to set it up experimen- 
tally. A suitable Hamiltonian must not only serve the 
given evolution mathematically, it should also be realiz- 
able physically. In particular, the problem becomes very 
difficult for nonunitary evolutions of open systems albeit 
its solvability for unitary evolution of closed systems. 

The purpose of this paper is to put forward a general 
approach for finding the appropriate Hamiltonian that 
determines a given evolution of the systems. We restrict 
our discussion to two level system (the qubit) which is 
generally the most prevalent system used in quantum 
computation and quantum information. We will pro- 
vide a general formalism of the Hamiltonians for realiz- 
ing an arbitrary prescribed evolution of the qubit system. 
The Hamiltonians are not unique. Our formalism gives 
a set of the equivalent Hamiltonians. Both unitary and 
nonunitary evolutions are investigated respectively. 



HAMILTONIANS FOR REALIZING UNITARY 
EVOLUTION OF A QUBIT 

For a qubit system, the density matrix can be ex- 
pressed as 



1 



p(t) = -(l + r.(T), 



(3) 
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where r, the Bloch vector, is a time-dependent vector 
function with |r| < 1 and a are the Pauli matrices. When 
we refer to a given evolution, it means that p(t), as a 
matrix function of t, is explicitly given. In other words, 
the movement of the Bloch vector r in Bloch sphere is 



known. The most general form of p{t) can be written as 
1 /l 



Pit) = 



r cos o r sm a e 
r sin 9 e"^ 1 — r cos ( 



(4) 



where r ~ r{t), 9 = 9{t) and (p — (j){t) are arbitrary time- 
dependent real functions with 0<r<l, 0<6'<7r, 0< 
(/) < 27r. At i = 0, we denote r(0), 61(0) and 0(0) as ro, 6*0 
and 4>o respectively, and the general initial state reads 



1 / 1 + ro cos 9o ro sin 6*0 e-"^° 
^^ ' 2 I ro sin 9o e"l"> 1 - ro cos 9o 



(5) 



Our aim is to find the general form of Hamiltonians that 
realize the evolution defined by Eq. 10} with the initial 
state Q. For unitary evolution, r = ro is constant and 
there exist unitary operators U{t) satisfying 



pit) = U{t)p{Q)U{t)l 



(6) 



If we can find the operators U{t), the Hamiltonians H{t) 
are easy to be obtained as 



H{t) = iU{t)U{t)K 



(7) 



Hence the problem reduces to the question of how to 
find the unitary operators U{t) satisfying Eq. ©. As 
mentioned before, the corresponding operators are not 
unique. There are infinitely many unitary operators U{t) 
that can realize the same evolution, all of which are equiv- 
alent for p{t). We can find all the different but equivalent 
Hamiltonians by finding all the U{t). To this end, we first 
construct one of the unitary operators. We find that the 
following operator 



U{t) = 



■ <H-<h) 



- sm 
cos 



■ .^ + <^o 



2 
0-00 



(8) 



satisfies Eq. (jHJ for an arbitrary evolution defined by Eqs. 
Q and 10. With the operator U{t), all the equivalent 
unitary operators can be constructed by 



U{t) = Lf(t)V{t), 



(9) 



where V{t) e U{2) and satisfies the following commuta- 
tion relation, 



[V{t) , p(0)] = 0. 



(10) 



To provide an explicit form for matrix V{t), we rewrite 
the initial density matrix p(0) in the form of the orthog- 
onal decompositions. 



p(0) = J2 ^kPkiO), 



(11) 



k=l,2 



where wi — {1 + ro)/2, W2 = {I — ro)/2 and 



1/1 + 4 cos 9o 6k sin 9o e-'^« \ ,^ „. 

P''^^^-2\5ksin9oe^^o i„4cos0o j' ^^^^ 
with (5i = 1, ^2 = — 1- The above orthogonal decompo- 
sition is unique if ro 7^ (we need not consider ro = 
because in this case the Bloch vector does not move at all 
during the unitary evolution). Thus, the unitary matrix 
Vit) satisfying Eq. l(Tn|) can be written as 



V{t) = J2 e'"Vfe(0), 

fe=l,2 



(13) 



where ak — ak{t) {k — 1, 2) are arbitrary real parame- 
ters with ai \t=Q— Oi2 |t=o= 0. From Eqs. @ and H13|l . 
the general form of the unitary operators V(i) is obtained 
as 



C/(t) = C/(f)^e'"V.(0). 



(14) 



fc=i 



The corresponding general form of the Hamiltonians Hit) 

is 



ff(t) = z?7(i);7(i)t-^d,[/(t)pfe(0)i7(t)t 

k=\ 
2 

= H(t)-Y,dikPk{t), (15) 



k=\ 



where iJ(i) = iU{t)U{ty and pfc(i) = f7(t)pfc(0)L/(i)t. 
More explicitly, the Hamiltonians read 



H{t) 



~ di(l + cos( 
{i9 — ai sin 6 



— a2{l — cost 



a2 sm9)e 



i(jj 



{-i9 - ai sin 9 + 02 sin 6')e-*'* 
—0 — ai{l — cos9) — 0:2(1 + cos( 



(16) 



r 



The formula expressed by Eq. (|16|l provides the recipe scribed evolution defined by Eq. 10} and |(S}. Since the 
for finding all the Hamiltonians that can realize the pre- 



evolution given by Eqs. Q and (O is arbitrary, the 
Haniiltonians given in the above equation are the most 
general ones. When a prescribed evolution p{t) is given, 
i.e., r = rg is a known constant, and 0{t), <p{t) are two 
known time dependent functions, one can easily write 
down the corresponding Hamiltonians by directly substi- 
tuting the variables into Eq. H16() . Noting that there are 
two time-dependent parameters ak{t) (fc = 1, 2), which 
can be arbitrary real functions of t with afc(O) = 0, one 
can establish the suitable Hamiltonian for any special 
purpose by properly choosing the forms of ak{t). 



variable r. Suppose the ancillary qubit is initially in the 
state |0)(0|, the initial state of the combined system can 
be written as q{Q) = p{Q) ® |0)(0|, where p(0) is an ar- 
bitrary given initial state of the open system. Then, at 
time t, the density matrix Q{t) of the combined system 
reads 



Q{t) = Uab{t)piS))®\Q){0\Uab(t)^- 



(17) 



The reduced density matrix of the open system is given 
by tracing out the state of the second qubit, i.e. 



HAMILTONIANS FOR REALIZING 
NONUNITARY EVOLUTION OF A QUBIT 

We now turn to the nonunitary evolution of an open 
qubit system. For nonunitary evolution, the path traced 
out by the Bloch vector is a curve inside the Bloch sphere 
with variable r. In this case, Eq. (Q is no longer valid. 
To realize a prescribed evolution of the open system, we 
need to consider a larger closed system, comprising the 
open system and an ancillary system which can also be 
regarded as the environment of the open system. For sim- 
plicity, we take the ancillary system to be a qubit. We 
denote the qubit under prepared as a and the ancillary 
qubit as h, respectively. The two qubits constitute a com- 
bined system Sab- The combined system is a closed one, 
and its evolution is described by unitary operator Uab{t). 
The most general evolution of an open qubit system can 
still be described in terms of expressions Q and Q with 



p{t) = tYbg{t) = ti-b{Uab{t)p{Q) <E> |0) {0\Uab{t)^}. (18) 

The problem of preparing an open qubit system for a 
given evolution reduces to that of finding the unitary op- 
erator Uab{t) such that the reduced matrix p{t) remains 
the same as that assigned by Eqs. Q and jSJ). To this 
end, we rewrite Eq. H18|l by using the Kraus representa- 
tion 

p{t) = Y, M^p{0)M^\ 
n=os 

J2mIm,,=I, (19) 

A' 
where the Kraus operators A/^ = {p\Uab\0) , and \p) [p = 
0, 1) are the orthogonal bases of the ancillary qubit. If 
we know the Kraus operators, we may inversely deduce 
the unitary operators Uab{t). One readily verifies that 
the following operators Af^ satisfy Eq. H19|) . 



Mo - 



Ml 



ro 



1 + r-o 



■^ sin I cos %e* 

l+ro 2 2 



1+ro 



cos ^ cos ■ 



cos § COS fe"l"> 



So_pi4> 



^^ sin I sin ^e"*"^ 

l+ro 2 2 



—^ COS ^ sm -^ 

\-ro z z 



(20) 



r 



The choice of the Kraus operators is not unique. There 
are infinitely many (equivalent) operators M^ satisfy- 
ing Eq. (|19|l for given p{t) and p(0). M^ is only 
one special choice. Correspondingly, there are also in- 
finitely many choices of Uab{t) satisfying Eq. 118|) . so do 
H{t). With the known Kraus operators M^{t), we can 
first construct one of the operator Uab{t), and denote 
it as Uab{t). In fact, Uab{i) is a 4 x 4 unitary matrix, 
where the elements in the first and the third columns are 



completely determined by the Kraus operators Af^ with 
f>ii = (A/o)ii, U21 = (Mi)ii, (731 = (A~fo)2i, Uai = 

(Afl)21, C/13 = {Mo)l2: U23 = (A^/l)l2, U33 = {Mo)22, 

and U43 ~ (Afi)22 while the elements in the other two 
columns are yet to be determined. The second and 
the fourth columns' elements can be chosen by ensur- 
ing Uab{t) to be unitary. Without loss of generality, we 
choose tJab{i) as 



Uab{t) 



/-Cs- 
V 



r^ScC 
r+CcC 



i{(f>o-<ti) 



Sse'f + r^Cce' 



r+Sce 



i{<i>+M 



-SnC 



CcC 



r+Sse 

''■{'t'o-tt') 



r-Cs 



'+<-s 



r-SsC 



i(j) 



SrC- 



-i(p 



~r-Sse '"^ 


r+Sce '^y'^^^" 


r+Cs 


—SgC'^'^ + r^CcC 


-M +r_cs 


-r+Cce~*'^« 



r+SsC 



i(p 



r-ScE 



i{4>-4>o) 



(21) 



where Sg = sin | sin ^ , Sc = sin | cos ^ , c^ 
cos I sin ^, Cc = cos I cos ^ and r+ = w-ji^, r_ 



l-r 



Expression (|21|) is only one of the many unitary opera- 
tors Uab{t) satisfying Eq. I|18() . To construct other equiv- 
alent unitary operators, we need to find other equivalent 
Kraus operators Af^. In terms of the known operators 
given by Eq. (|2()|l . the equivalent Kraus operators can be 
written as 



M^ = 



W^,M,V{t), 



(22) 



where W^^ are the elements of time dependent matrix 
W{t) e SU{2), and V{t) € U{2) satisfies Eq. ^. 
The explicit expression of V{t) is given by Eq. H13() if 
rg ^ and it can be any arbitrary 2x2 unitary ma- 
trix if rg = 0. Eq. H22|l is the general expression of the 
Kraus operators realizing the prescribed evolution of the 
open system. With the general Kraus operators M^ , 
we can construct the general expression of unitary oper- 
ators Uab{t)- By using Eq. (|22|l . the equivalent set of the 
unitary operators can be described by 



Uabit) = il(g> W{t))Uabit)iV{t) (g> /), 



(23) 



Substituting Eq. H23I) into Eq. iQ, one obtains the 
Hamiltonians 



Habit) - iUab{t)Uab{t)l 



(24) 



For the prescribed evolution of the open system, there are 
many possible choices of Hamiltonians, one of which is 

given by Hab{t) = iUab[t)Uab{i)'^ while the others are ob- 
tained by Eq. (|^ . For a given prescribed evolution p(i), 
we can immediately write down the Hamiltonians by sub- 
stituting the corresponding r — r{t), = 9{t), (p — (j){t) 
into Eqs. ^, ^ and (EUl. AU the Hamiltonians fulfill 
the prescribed evolution of the open system. As the ma- 
trices W G SU{2) and V{t) = J2 e'"'=/5fe(0) if tq ^ 0, 

fe=l,2 

or V(t) G U{2) if ro = 0, there are five or seven pending 
parameters to serve any other special requisitions. One 
can obtain any required Hamiltonian by choosing prop- 
erly the pending parameters. 



EXAMPLE FOR APPLICATIONS 

In the previous sections, we develop a general proce- 
dure of Hamiltonians for realizing a prescribed evolution 
of a qubit system, be it closed or open. In spite of the 
complicated details, it is still a useful device due to its 
generality. As an illustration of its applications, we con- 
sider the measurement of geometric phase as an example. 

As we know, when a system evolves along a given path 
in the project Hilbert space, the system will acquire the 
dynamic phase as well as the geometric phase. Geometric 
phase can be regarded as the total phase minus dynamic 
phase. For the mixed state of a qubit system, if the 
evolution is defined by Eqs. Q) and ISJ with r = tq, the 
geometric phase acquired by the system is 3, 4] 

where wi == (1 + ro)/2, W2 = {I — ro)/2, and Pfe(O) are 
defined by Eq. 1)12(1 . Mathematically, we may use any 
one of the unitary operators U{t) given by Eq. p4|l to 
evaluate the geometric phase. All of them result in the 
same value of j{t). However, the experiment situation 
for phase measurement is quite different. The phase mea- 
sured in the experiment is usually the total phase of the 
system, not the geometric phase. The different choices 
of U{t) will lead to different values of the measurement. 
If one wishes to measure the geometric phase, one must 
choose a special unitary operator satisfying the parallel 
transport conditions 



tr[pfe(0)C/(f)t{7(t)]=0, (fc-1,2) 



(26) 



The special unitary operator will render the dynamical 
phase zero, hence the total phase equals the geometric 
phase. Substituting Eqs. ifT^ and (fTljl into Eq. if^ . one 
will find that the parallel transport conditions constrain 
the functions ai{t) and 012 (i) to be [5j 



Qfl = 



-a2 = - 



cos 



(27) 



Substituting Eq. lf?7|) into Eq. (fTB|l . the Hamiltonian is 
obtained as 



Hit) = iU{t)U{ty = 2 ^ (^g 



ai4> 



')sm6cos6)e 



(28) 



r 



with which the measurement for the total phase results in 
the geometric phase. Let us be more specific. Suppose a 
qubit system, say an electron, with initial state p(0) (tq ^ 
1, 00 = 0) given by Eq. jSJ, is expected to evolve along a 
path with 9 = 9o, 4) = tot , t € [0, ^]. The evolution is 
a closed circle inside Bloch sphere. By using expression 
(|28|l . the special Hamiltonian can be immediately written 
down as 



Hit) 



Lu sin 9o 



smt/o 
— cos 00 e" 



cos ooe 

— sin 9o 



(29) 



The above Hamiltonian ensures the given state to evolve 
along the prescribed path with zero dynamic phase. The 
measurement result of the geometric phase should be 
7 = — tan~-'^{ro tan[7r(l — cos6'o)]}. This Hamiltonian 
can be experimentally realized in physics by applying a 
time-dependent magnetic field, B(t) = (_Bi, B2, B^) — 
a;sin0o(~ cos^o coscjt, — cos^q sintji, sin^o); to the 
electron system, and Hit) — f7-B(t)/2. As the geometric 
phase has been used as quantum gate, the above discus- 
sion may be useful for quantum computation. As well- 
known, if an evolution causes nonzero dynamic phase, 
one must eliminate it by some means[2(-|a|- However, this 
cancellation is not easy. In our case, given the generic ex- 
pression for the Hamiltonian, one can try to select those 
Hamiltonians for which the dynamic phase is zero. Recall 
that Eqs. H28() and (|29|) give no dynamical phase for any 
path of the state, be it closed or open. It is also inter- 
esting to consider non-cyclic evolutions instead of cyclic 
closed circuit. In this context, we would like to mention 
a recent paper 9], where it is argued that some systems 
may obey the same fault-tolerance properties like those 
of the geometric phase also for nonvanishing dynamical 
phase. 

The above example provides a concrete application of 
our formalism in the form of the geometric phase of mixed 
state under unitary evolution. Similarly, we can also con- 
sider examples of mixed state under nonunitary evolu- 
tions. For instance, suppose we need to prepare a qubit 
system whose Bloch vector r is expected to evolve along 
a prescribed path inside the Bloch sphere, say, an ellipse 
with r\ -I- 4r^ = 1, r^ = 0. We can realize this evolution 
through the combined system of the qubit and an ancil- 
lary qubit. The unitary operator Uabit) can be easily cal- 
culated by substituting r — (cos^ ujt + A sin^ ujt) '^ .9 — 
|, = uot, and ro = 1, 6*0 = f , 0o = into Eq. (EU- 
With Uabit), the general unitary operators Uab{t) and 
the general Hamiltonian Habit) can be obtained by Eq. 
(P^ and Eq. (|^ respectively. Each Hamiltonian Hit) 



obtained in this way satisfies the same evolution with 
the appropriate reduced matrix pit) along the given el- 
lipse. If there are other physical constraints, one can 
then choose a preferred Hamiltonian by determining the 
suitable parameters in the Hamiltonians. As the Hamil- 
tonians involve complex 4x4 matrices again, we will 
not describe the example in detail in the present paper. 
Instead, we prefer to show an explicit form of the Hamil- 
tonians for a simple example, the pure shrinking of the 
Bloch vector along the polar axis. In this case, by putting 
6* == 6*0 = 0, To = 1 in Eq. if^ . the Hamiltonian Habit) 
is obtained as 



Habit) = iUab{t)Uabit) 



t _ 



4vr 



. a^l^O) 



where a^, (Jy are Pauli matrixes, and r (0 < r < 1) is an 
arbitrary prescribed time dependent function with f < 0. 
Habit) will render the Bloch vector of the state shrink 
from the pole (1,0,0) towards the center of the Bloch 
sphere with the prescribed speed \f\. Other Hamiltonians 
equivalent to Habit) can be obtained by using Eqs. H23(l 
and if^ . 



DISCUSSIONS AND CONCLUSIONS 

The inverse problem for the evolution of a qubit is 
investigated in this paper. We proposed a general for- 
malism to establish the Hamiltonian for a qubit system 
that is required to evolve along a particular path. Both 
unitary and nonunitary evolutions of mixed states of a 
qubit system are discussed. 

For a closed qubit system, its evolution is unitary. The 
general form of Hamiltonians for realizing a prescribed 
evolution is explicitly given by Eq. (|16|l . One can directly 
write down the Hamiltonians by substituting 9 — 9 it) 
and 4> = 'Pit) into Eq. (|16|l . The explicit expression for 
the Hamiltonian contains two arbitrary time dependent 
parameters. Regardless of the values of the parameters, 
all Hamiltonians give rise to the same evolution. One 
may choose appropriate parameters to satisfy the physi- 
cal constraints. 

For an open qubit system, its evolution is generally 
nonunitary. However the system with its environment 
continues to evolve as a closed system. We therefore in- 
vestigate the open system by considering the closed sys- 
tem constructed by adding an ancillary qubit to the open 
system. The evolution of the open system is expressed 
by the reduced density matrix given by Eq. H18|) . If the 



desired evolution is as given by Eqs. Q and (0), the uni- 
tary operators Uab{t) and the Hamiltonians Hab{t) can 
be calculated by Eqs. H23|) and 1)24(1 . The expression of 
the Hamiltonians contain five or seven arbitrary parame- 
ters. All Hamiltonians of the given form cause the qubit 
to evolve along the desired path while they really act 
on the combined system. One may choose the suitable 
parameters for the Hamiltonians to meet any physical 
constraints, just like the unitary case. 

As an example for illustrating the applications of our 
result, we discuss the parallel transport conditions on 
geometric phase measurement. A general expression of 
the Hamiltonians satisfying the conditions is given by 
Eq. H28|l. With the kind of Hamiltonian, the value of 
the measurement on total phase results in the geometric 
phase. 

The study in present paper focuses on the qubit sys- 
tem. However, the approach can be generalized to higher 
dimensional systems. Some main formulae are similar. In 
fact, for unitary evolution of A^-dimension system, Eqs. 
(fO|l - (fT5| are stiU valid with fc = 1, 2, ..., A^. For nonuni- 
tary evolution of Af-dimension open system, Eqs. (|23|) 
and H24|) still hold. The difference for an A-dimension 
system from the qubit system appears in U{t) given by 
Eq. © or Uabit) given by Eq. ij^lj, which wiU be re- 
placed by some new matrices corresponding to an A^- 
dimensional system. 
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